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We present numerical results on calculations of energy spectra and wave functions of annular 
billiard. 


I. INTRODUCTION 


In recent decades a new direction of investigations 
started its active developing devoted to studies of the 
so-called quantum chaos. This term implies analysis of 
behavior of the quantum systems whose classical coun¬ 
terpart demonstrates chaotic dynamics. 

Among the systems demonstrating chaoticity, one of 
the simplest family is presented by plain billiards, which 
are two-dimensional single-particle potential-free sys¬ 
tems, limited in finite domain. They are defined by two- 
dimensional single-particle Hamiltonian 


H = 


2 , 2 
P x +P y 

2m 5 


for r £ fl; 


oo, for r ^ Q. 


The dynamics (both the classical and quantum) in such 
systems is defined by solely the geometry of the bound¬ 
ary. In the present paper the boundary is chosen to have 
the following form: motion takes place in the region lim¬ 
ited by two circles of different radii (see FigjT]). Such 
system is called the annular billiard. In the case 5 = 0 
(centers of the bounding circles coincide) the system has 
two independent integrals of motion: energy and angular 
momentum. As motion takes place in two-dimensional 
space, it is sufficient to make the dynamics purely reg¬ 
ular. For 5 > 0 (inner disk is shifted relative to the 
external one) the system looses rotational symmetry and 
the angular momentum is not already conserved, thus the 
system possesses mixed (chaotic and regular) dynamics. 
The annular billiard is a well-known system widely used 
in the present context m- Besides the convenience for 
analysis of fundamental properties of quantum chaos, the 
annular billiard represents interest for experimental and 
applied researches, in frames of which it models the mi¬ 
crowave cavities nig. 

A characteristic feature of chaotic system which man¬ 
ifests in quantum dynamics is the dynamical tunneling. 
This notion generalizes usual tunneling. The dynamical 
tunneling represents a transition between two quantum 
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FIG. 1: Boundary shape of the annular billiard. The param¬ 
eter 5 is considered to be positive. 


states which is forbidden for their classical counterparts 
for some reason (not obligatory because of static poten¬ 
tial barrier) [g. 

In 1993 year Bohigas and colleagues |J] were first to 
describe the dynamical tunneling in the annular billiard. 
There is a family of trajectories in this system which do 
not touch the inner disk. They are usually called the 
whispering gallery modes. As those trajectories do not 
feel the inner disk, they are regular and quasiperiodic. In 
his paper Bohigas demonstrates the transition between 
the two states, corresponding to two whispering gallery 
orbits. This phenomenon has a distinct peculiarity: the 
tunneling rate varies in many orders of magnitude with 
comparably small changes of system parameters. 

One more important feature of chaotic quantum sys¬ 
tems is the presence of level quasicrossings (level repul¬ 
sions, avoided crossings). The theorem proven by Wigner 
and von Neumann [7j states that two initially different 
values of observables, which correspond to states in the 
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FIG. 2: Example of the level repulsion. A typical dependence 
on a parameter of Hamiltonian is shown for two energy levels 
of the states of identical symmetry type. 


same symmetry type, can be made equal by variation of 
Hamiltonian parameters at most on a manifold of dimen¬ 
sion N — 2, where N is number of the parameters varied. 
Therefore variation of only single parameter cannot lead 
to intersection of energy values of two levels and thus it 
forms a hyperbola in vicinity of the nearing point (see 
Fig{2]). The quasicrossings are related to such properties 
of quantum chaos as changes in the level statistics and 
chaos-assisted tunneling [2]. It was one of the manifesta¬ 
tions of the latter effect which was described in the Bo- 
higas’ work. There remains an open question about the 
details of quantum-classical correspondence in the case of 
the mixed type of dynamics. It follows from qualitative 
considerations that the wave functions of the states in 
such systems corresponds to some invariant sets of clas¬ 
sical dynamics. However such assumptions does not hold 
near the quasicrossing points on the intermediate stage 
of the mode exchange. In the present work we make 
a detailed analysis of localization features for (quasi- 
distributions in the phase space of quantum states on 
the classical trajectories of the system. For that purpose 
we use the concept of the Husimi-Poincare functions. We 
consider a quasicrossing between a regular and a chaotic 
state. Besides that we explicitly demonstrate tunneling 
between them near the quasicrossing point. 


II. PROBLEM SETUP AND NUMERICAL 
METHOD DESCRIPTION 

A quantum analog of classical billiard represents a 
particle freely moving in a region fl bounded by elastic 
walls. It is formally expressed by the following stationary 
Schrodinger equation 


yZ 

Eip( r) = —— AV>(r), red, 

2m 

with the Dirichlet boundary condition 
4 >{r G dtt) = 0, 


(1) 

( 2 ) 


where dfl denotes boundary of the region Q. 

Non-trivial shapes of boundary of the domain H do not 
allow analytical solutions of the equation 0 - Therefore 
it is necessary to apply numerical methods. First for con¬ 
venience let us rewrite the equation (nj , substituting the 
energy by the wave number squared (below and further 
we use the system of units with 2 m = h 2 = 1): 

( k 2 + A)ip(r) =0, r G Cl. 

There are plenty of numerical methods for solution of 
Schrodinger equation, however not all of them can be ef¬ 
ficiently applied to quantum billiards. The main feature 
of such systems—absence of the potential—can be used 
for elaboration of efficient algorithms for the numerical 
solution. A widely used family of such algorithms is pre¬ 
sented by the boundary integral method. Its efficiency is 
based on utilization of the fact that the eigenfunction of 
the equation 0 is completely determined by the values 
of normal derivative on the boundary of the billiard. Let 
us introduce the normal derivative as a function of the 
natural length parameter s along the boundary: 

u(s) = r(s)), r(s) G dtt. 

The wave function in arbitrary point is expressed through 
the normal derivative on the boundary in terms of the 
Green function of the equation (flj) in the following way 
0 : 


ip{r) = I Gfc(r, r(s))u(s)ds. 
an 

Following the paper [5] we may transform Schro’odinger 
equation for our quantum billiard into an integral equa¬ 
tion by means of the regularized Green function G&( r, r') 
which solves the following defining equation: 

(A + k 2 )G( r, r') = <$(r — r') — d(r — r^), (3) 

where r and r' are in B U dB , and is the mirror image 
of r' with respect to the tangent at the closest lying point 
on the boundary, and thus if r' is sufficiently close to the 
boundary then r R is outside the billiard B. The solution 
can easily be found in terms of the free propagator (the 
free particle Green function on the full Euclidean plane) 

G 0 (r,r / ) = (4) 

where H^ 1 ' 1 is the zero order Hankel function of the first 
kind m, namely 

G(r, r') = G 0 (r, r') — G 0 (r, r^), (5) 

such that now G(r,r') is zero by construction for any r' 
on the boundary. 
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With the use of the normal derivative the equation (jTJ) 
with the boundary condition Q can be rewritten in form 
of the integral equation 

u{s) = -2 -^Gk(r(s),r(s'))u(s')ds'. (6) 

on 


For the normal derivative of the Green’s function one 
obtains 


d_ 

dn 


G k 


ik n(s) • (r(s) — r(s')) 
4 |r(s)-r(s')| 


x 7J^(&;|r(s) 


r(OD- 


sv 



— SV1 

— SV2 


Therefore the two-dimensional problem is reduced to the 
one-dimensional one which is much easier to solve numer¬ 
ically. For that purpose we present the boundary integral 
in the equation © in form of finite sum. We thus obtain 
a finite system of uniform linear equations 

y =o, (7) 

k 

where where is a dense, complex non-Hermitean ma¬ 
trix [5]. Solution of the equation ([7]) requires to find those 
values of the wave number k which turn the system de¬ 
terminant det(A) to zero. It is the values of k that give 
the required stationary energy levels. 

A difficulty arises here: there are no real values of k 
which make the obtained determinant equal to zero ex¬ 
actly, because it is a complex function of a complex ar¬ 
gument. Due to the discretization of the integral © the 
determinant det(A) will not become zero but only close 
to zero m- Besides that, calculation of the wave func¬ 
tions requires to find the normal derivative u(s) in the 
discretizing points of the boundary. To get rid of this 
problems one should in advance to make the singular de¬ 
composition of the matrix A 

A = USV\ 


FIG. 3: Dependencies of first and second singular values on 
the wave number in symmetric annular billiard (r = 0.2). One 
can see that near the value k = 77 both the singular values 
have a sharp minimum, which gives evidence for presence of 
a multiplet. 
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FIG. 4: The first singular value calculated with improved 
accuracy is shown here in the vicinity of minimum of the 
second singular value (see Fig.3). It can be seen that one can 
reveal two values of the wave number k by decreasing of the 
calculation step and increasing of the accuracy. 


where U and V are unitary matrices, and A is a diago¬ 
nal real matrix composed of singular values of matrix A. 
Those singular values are defined as eigenvalues of the 
self-adjoint operator AA*. Minima of the first (minimal) 
singular value as function of the wave number correspond 
to the sought eigenvalues, and the corresponding columns 
of the matrix U give values of the normal derivative. Be¬ 
sides that, such approach allows to detect the energy val¬ 
ues which lie near the quasidegenerate states. If both the 
first and second singular values have sharp minima near 
the obtained eigenvalue then there are in fact two states 
in close vicinity of the energy eigenvalue. 

Examples of calculation of the singular values are pre¬ 
sented on Figures [3] and [4] One can easily see from this 
example that usage of the singular decomposition allows 
to detect the quasicrossings and find all the quasidegener¬ 
ate levels with considerable saving of the computational 
efforts. First of all one should find all minimum values 
of the first singular value. In the case of detection of 


minimum for the second singular value one repeatedly 
performs search for required values of the wave number 
with ever decreasing step and higher discretization degree 
of the matrix A , aimed to improve the accuracy. To ver¬ 
ify the obtained results we use the level statistics. Let 
us introduce the integral level density—the state num¬ 
ber staircase function, which counts the number of levels 
with energy not greater than the given one 

N(E) = {n\Ei< n < E, E n+1 > E}. (8) 

For billiards one can use the generalized Weyl formula 
which gives the averaged smooth component of the func¬ 
tion ® [II] : 

N ^ E )=^-AVe+C. (9) 

Here we dropped the terms of higher orders of smallness, 
S is the interior area of the domain f 1, L = Ljj — Ljy, 









4 


iV c (E)-iV w (E) 




FIG. 5: The fluctuating part of the cumulative distribution 
function for symmetric (upper) and antisymmetric (lower) 
states obtained numerically. 


Ld/n are the parts of the boundary with the Dirich- 
let and Neumann conditions respectively. The desym- 
metrization we use here allows to consider separately 
the level statistics for states of different symmetry. At 
all that one should take into account that the part of 
the boundary along the symmetry line of the billiard 
has Neumann boundary conditions for symmetric states 
and Dirichlet boundary conditions for the antisymmetric 
ones. Subtracting the smoothed part © from the numer¬ 
ically obtained integral level density N C (E), one obtains 
the fluctuating part. In the case when a level is missing 
in the numerical results it can be immediately detected 
by strong deviation from zero of the average value of 
the fluctuation component. Fig[5]shows the obtained de¬ 
pendencies, which give evidence for completeness of the 
calculated spectrum. 


III. QUASICROSSINGS OF QUANTUM STATES 
AND DYNAMICAL TUNNELING 

Consider an initially regular annular billiard with S = 
0. Let us start gradually shift the inner disk, i.e. to in¬ 
crease S. The tori corresponding to orbits near the inner 
disk will be destroyed and distorted. The quantum states 
corresponding to them will get all their properties con- 



FIG. 6: Dependence on the shift parameter S for two en¬ 
ergy levels of antisymmetric states. It is an example of level 
repulsion in the annular billiard. 


siderably changed, in particular their energy significantly 
changes with variation of 5. However as Lazutkin have 
shown P2], in the case of smooth billiard boundary[TT] 
there are non-destroyed tori remaining in the system for 
any set of parameters. In the considered case of the an¬ 
nular billiard there are whispering gallery orbits up to 
S = 1 — r. As the whispering gallery tori in classical dy¬ 
namics absolutely do not feel the variation of the 6 (up to 
certain its value), it is then natural to assume that the 
semiclassical quantum states corresponding to the tori 
will be extremely weakly sensitive to the variations of <5 
too. 

Values of some energy levels get closer at variation of 
5. If these levels have the same symmetry then after 
certain nearing they start to repel and form a quasicross¬ 
ing. While passing the quasicrossing region, these states 
exchange the so-called adiabatic properties. We are in¬ 
terested in peculiarities of this exchange in the case of 
quasicrossing between chaotic and regular levels. The 
simplest task is to detect such mode exchange between 
the states which correspond to classical orbits of whis¬ 
pering gallery and chaotic ones. This is due to the fact 
that, as was mentioned before, the energy levels of the 
former weakly depend on the shift parameter, while the 
letter depend much stronger. 

Figure [6] presents the dependence of energies on pa¬ 
rameter S for two antisymmetric states. The very mode 
exchange takes place in very narrow range of energy and 
parameter space (of order 10 -4 ). Thus numerical results 
with lower accuracy may look as real level crossing. Near 
the quasicrossing the energy level lines have characteristic 
hyperbolic shape. As can be seen, one of the levels has 
inherently practically constant energy value, the other 
considerably shifts with growth of 6. One can thus as¬ 
sume that the former state corresponds to the whispering 
gallery and the latter—to the state related to chaotic sea 
of classical dynamics. 

Figure [7] shows wave functions of the above described 
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FIG. 7: The same as on Fig.6 with addition of images of wave 
functions for some states (the corresponding energy values are 
shown by arrows). 


states. Let us first consider the wave functions of the 
states in relative distance from the quasicrossing. First 
of all it should be noted that the modes under investiga¬ 
tion completely exchange the wave functions. Distinction 
of the wave functions to the left and to the right of the 
quasicrossing is extremely weak and it is related to vari¬ 
ation of the parameter S. The state corresponding to 
the horizontal line has regular shape of the wave func¬ 
tion as expected. This wave function is localized near 
the outer disk, it strongly decreases near the inner one, 
which shows its connection to the whispering gallery or¬ 
bits. The other wave function, also as expected, is clearly 
chaotic, what is supported for example by irregularity of 
its nodal lines. More rigorous correspondence of wave 
functions of these states to the classical structures will 
get more detailed investigation in the next section. Wave 
functions near the quasicrossing point were calculated for 
parameter value which provides the minimum of the en¬ 
ergy splitting. It can be seen that the states are close to 
the original ones but they are mixed in some sense. The 
lower medium state on the Figure[7]is clearly close to the 
regular one, which is far from the quasicrossing point. 
Nevertheless its wave function is substantially distorted. 
There are traces of the chaotic wave function clearly vis¬ 
ible in the central part. 

It is natural to assume that during the mode exchange 
the wave functions represent linear combinations of the 
original ones. However this notion cannot be understood 
too literally here, because variation of the parameter <5 
affects the configuration space itself. But smallness of the 
parameter changes allows to verify such a hypothesis. Let 
us construct a linear combination of the wave functions 
near the quasicrossing point. The coefficients should be 
fitted in order to make the obtained combination simi¬ 
lar to one of the original wave function. For example, it 
can be seen from the Figj7]that the wave function corre¬ 
sponding to the whispering gallery mode has a nodal line 



FIG. 8: Consecutive stages of wave function evolution. 


along the vertical diagonal of the outer boundary. Mini¬ 
mization of the absolute value of the considered combina¬ 
tion on this line gives a combination which is practically 
identical to the original wave function[T8]. The state or¬ 
thogonal to the latter in turn recovers the original chaotic 
one. 

Having obtained the chaotic and regular combinations 
of the wave functions, one can observe the dynamical 
tunneling between these states. In the billiard with the 
parameter 5 close to the quasicrossing point we prepared 
the regular wave packet described above. Figure [8] shows 
the time evolution of this ave packet. As can be seen, 
the regular state (partially) tunnels into the chaotic one. 
It is initially manifested as the wave function, initially 
concentrated on the boundary, (partiallyO spreads over 
whole configuration space. Distraction of the initial ro¬ 
tational symmetry takes place. The nodal lines loose 
regularity, forming new structures and distorting the re¬ 
maining ones. After that the states returns to the original 
one. 


IV. PECULIARITIES OF 
QUANTUM-CLASSICAL CORRESPONDENCE 

We need to relate the quantum states to the classi¬ 
cal dynamics. Different quasi-distributions of quantum 
states in phase space represent a convenient tool for this 
purpose. It is impossible to construct actual distribution 
because of a number of restrictions imposed by the un¬ 
certainty principle. Therefore one have to use different 
variations in some sense similar to the distribution func¬ 
tion. In the present work we shall use the concept of the 
Husimi functions. It satisfies all the criteria required for 
a distribution functions except the axiom of additivity of 
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probabilities in terms of Kolmogorov axiomatic. 

The Husimi function is defined (up to a coefficient 
independent on phase variables) as the absolute value 
squared of the folding of the wave function with the co¬ 
herent state JH]. However it is problematic to work with 
complete Husimi functions because the visualization in 
four-dimensional space presents a difficulty. 

To solve this problem let us address the classical dy¬ 
namics. Absolute value of momentum of a particle, as 
well as its energy, does not affect shape of trajectory of 
its motion in a billiard. Only direction of the motion 
and position of the particle matters. Moreover motion of 
the particle is completely determined by position of its 
last collision with the boundary and value of the incident 
angle. Therefore motion of the particle can be reduced 
to discrete mapping in coordinate space defined by the 
pair (position of the collision with the boundary, cosine 
of the incident angle). This pair is called the Birkghoff 
coordinates. Such mapping represents in fact a kind of 
Poincare map. 

As was mention before, the quantum dynamics in 
billiards can be completely reduced to that of normal 
derivative of the wave function on the boundary. This 
fact shows an analogy with classical dynamics. Let us 
construct an analogue for Husimi function for the nor¬ 
mal derivative on the boundary—the so-called Husimi- 
Poincare function. For that purpose one has to make 
further improvements. First of all it is necessary to pro¬ 
vide periodicity of the function on the boundary, i.e. take 
into account its closure. For that purpose we introduce 
an one-dimensional coherent state which is periodic on 
the outer boundary m 


ck(q,p\s) = - 

7r 


E 


1/4 oo 

o ik\p(s— g+27rm) + ^ ( 5 — q-\-2nm) 2 ] 

(10) 

where k is the wave number of the state, which the 
Husimi-Poincare functions is calculated for; s is the 
natural length parameter along the boundary; (q,p) £ 
[0,2tt) x R are arguments of the function. 

Calculation of folding of the two functions on the 
boundary is quite arbitrary as the metrics of the two- 
dimensional domain does not induce any metrics on the 
boundary. Let us define the scalar product on the bound¬ 
ary in the simplest way: as an integral with unit kernel. 
Taking projection of the normal derivative of the wave 
function for a given state u(s) with the coherent state 
(10) one obtains the Husimi-Poincare function 


^ = 2 Tk 


2tt 


ct(q,p; s)u(s)ds 


The coordinates (q,p) represent direct analogy t the 
Birkghoff coordinates. Having construct the classi¬ 
cal Poincare map and compared it with the calculated 
Husimi-Poincare functions, one can determine which 
classical object correspond to certain quantum state. 



FIG. 9: Poincare surface of section (the stroboscopic map) 
for the annular billiard. 



FIG. 10: The Husimi-Poincare distribution functions for 
two states (the left ones on Fig|rj) mapped over the classical 
Poincare surface of section. 


A typical Poincare map for the parameters correspond¬ 
ing to the quasicrossing under investigation is presentad 
on Fig. [9] Upper and lower domains filled with horizon¬ 
tal lines are the whispering gallery orbits. The central 
part corresponds to the chaotic sea with several regular¬ 
ity islands. Now we can relate the classical and quantum 
mechanics. Figure 10 shows Husimi-Poincare functions 
for regular and chaotic states far from the quasicross¬ 
ing. The Poincare surface of section is overwritten on 
them. As was expected, one of the states is localized [20] 
on the whispering gallery orbits, and the second —in the 
chaotic sea. Figure m gives an analogous picture for 
one of the quasicrossing levels. It is easy to see that the 
Husimi-Poincare function partially lies in the chaotic sea 
and partially- in the whispering gallery orbits region. It 
should be noted that in the chaotic part the function is 
practically zero on several stability islands and in general 
almost does not touch them. Transition of the function 
into the whispering gallery region occurs in two steps. 
The first to penetrate into this regular region are those 
parts from the chaotic sea which were initially close to it. 
After that the regular regions start to widen until they 
transform to solid line similar to that shown on the Fig¬ 
ure [l0|a.. At last this solid line is joined by the remnants 
localized in the chaotic sea. 
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FIG. 11: The Husimi-Poincare distribution function for a 
state near the level repulsion point (the upper middle state on 
Fig0 mapped over the classical Poincare surface of section. 
Separate solid lines mark the boundaries of the whispering 
gallery orbits region and regular islands in the chaotic sea. 


V. THE CONCLUSIONS 


We used a numerical method to obtain energy levels 
and wave functions of the quantum states in the annu¬ 
lar billiard of the region of spectrum where chaotic ef¬ 
fects show up. Variation of the inner disk shifting pa¬ 
rameter in the billiard give examples of repulsion be¬ 
tween chaotic and regular levels. Comparison of quan¬ 
tum functions of the Husimi-Poincare quasidistributions 
and classical Poincare surfaces of section allowed to per¬ 
form detailed analysis of the correspondence between the 
obtained quantum states and classical orbits. For the re¬ 
pelling states we obtained the time evolution of the dy¬ 
namical tunneling. 


The demonstrated features of behavior of the quan¬ 
tum states near the quasicrossing point present interest 
for future studies. The non-trivial character of the local¬ 
ization of quantum states on classical objects and transi¬ 
tions of wave functions during the mode exchange shows 
that there are many problems along these directions still 
requiring proper solutions. 
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